Minimal resonances in annular non-Euclidean strips 
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Differential growth processes play a prominent role in shaping leaves and biological tissues. Using 
both analytical and numerical calculations, we consider the shapes of closed, elastic strips which 
have been subjected to an inhomogeneous pattern of swelling. The stretching and bending energies 
of a closed strip are frustrated by compatibility constraints between the curvatures and metric of 
the strip. To analyze this frustration, we study the class of "conical' closed strips with a prescribed 
metric tensor on their center line. The resulting strip shapes can be classified according to their 
number of wrinkles and the prescribed pattern of swelling. We use this class of strips as a variational 
ansatz to obtain the minimal energy shapes of closed strips and find excellent agreement with the 
results of a numerical bead-spring model. We derive and test a surprising resonance condition for 
strips with minimal bending energy along the strip center line to exist. 
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The elastic buckling of thin sheets plays an important 
role in the shaping of biological tissues, especially the 
leaves of many plants [i|-l5|. In this context, the buckling 
results from the addition of extra material in the sheet, as 
when the sheet is subjected to local growth [6]. Recent 
experiments have also demonstrated that such growth 
processes can be mimicked by synthetic systems jS-Q- 
There have been several different theoretical formulations 
used to predict shapes of buckled objects which, broadly, 
fall into two classes. The first approach is to incorpo- 
rate swelling b y rn odifying the Foppl-von Karman equa- 
tions 0, 0, [i3, im • This is especially suited to studying 
the stability of nearly flat sheets but its applicability for 
studying the post-buckling behavior of swelled sheets is 
not obvious. A second approach defines a metric, pre- 
scribed by the local swelling of the sheet, from which 
an "ideal," strain- free shape can be determined 12-18]. 
This approach is best suited for thin sheets, where one 
can expect the sheets to be nearly strain-free. 

In this paper, we will continue an approach, begun in 
Ref. 18], to study narrow strips that have been sub- 
jected to an inhomogeneous pattern of swelling. We will 
focus on the specific case of closed, unknotted strips such 
as one that may be cut from the edge of a circular disk. 
Based on the results of a direct numerical minimization of 
a model of a swelled strip, we will focus on a class of strips 
for which the normal vector on the center line always lies 
tangent to a fixed cone: a "conical" strip. Within this 
class of conical strips, we are able to find profiles, depend- 
ing only on the metric tensor along the center line, that 
agree favorably with numerical minimization of a bead- 
spring model. This analysis extends the results of Ref. 
|18| to the entire class of conical strips, and corroborates 
those results numerically. 

For very thin, open strips, it can be shown that near- 
isometric embeddings exist both with vanishing mean 
curvature and mean curvature gradients along the center 
line [2g|. Nevertheless, we derive a "resonance" condi- 
tion, written in terms of the metric and the number of 



wrinkles in the embedding, for which a closed conical 
strip can have vanishing mean curvature along its center 
line. Our results indicate that the additional topological 
condition that the strip be closed results in nonvanishing 
mean curvature ~ thus, topology frustrates the bending 
energy. This prediction is in agreement with our numer- 
ical results. 

Throughout this paper, we will specialize to the case 
that the Gaussian curvature is negative. As a moti- 
vational example, consider the experiment described in 
Klein et al Q, in which a thin film of NIPA polymer 
gel is endowed with a controlled, axisymmetric monomer 
concentration that varies as a function of distance from 
the disk center. When the temperature is increased, the 
disks shrink to a configuration in which if < is approx- 
imately constant. 

It is fortuitous that there are families of smooth sur- 
faces that realize shapes of constant K < |19l - [2l| . 
and, naively, one might expect a sufficiently thin sheet to 
adopt one of these known shapes. The reality, however, 
is that the experimental sheets seem to adopt extremely 
wrinkled morphologies that are quite different from the 
known smooth embeddings. Though our results are not 
directly applicable to these experiments because they are 
for narrow strips, they do shed some light on the nature 
of the frustration. In particular, we find numerically that 
the strips allow localized stretching in order to lower their 
bending energy, despite the ready existence of embed- 
dings that admit no stretching at all. This detail will be 
the subject of a future publication |22| . 

In section HI we frame the problem of inhomogeneous 
swelling in terms of geometry and isometric embeddings. 
In section |lll we specialize to closed strips with inho- 
mogeneous patterns of swelling and introduce the class 
of "conical strips" . In section IIIIl we describe our di- 
rect numerical minimizations and compare with analyt- 
ical predictions. Finally, in section IIVI we discuss our 
results. 
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FIG. 1. (color online) Geodesic coordinates can be construct- 
ing starting from a curve (thick line) parametrized by f, at 
which u — 0. The geodesies perpendicular to the u = curve 
are parameterized by their arc length u (thin lines). The 
curves perpendicular to the geodesies (dashed lines) consti- 
tute a coordinate system in which the metric takes the form 
ofEq. (0. 
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u^dv^ , where u, v are the radius and polar angle in polar 
coordinates, respectively. 

Coordinate systems of these types always exist locally 
and can be found by the construction illustrated in Fig. [T] 
(also see page 286 of Ref. [2J|). First, choose an arbitrary 
curve on the surface to lie along u = 0, and let v be the 
parameter along that curve. For instance, this curve will 
be the center line of the strips that we consider. Curves 
of constant v are defined to be geodesies perpendicular 
to the u — curve, parametrized by their arc length u. 
Curves of constant u ^ are thus curves of fixed distance 
u from the u = curve along that family of geodesies. 
Since the resulting curves are mutually orthogonal and 
u is arc length, the metric must have the form of Eq. 
P]). By Gauss's theorema egregium, the Gaussian cur- 
vature can be written in terms of the metric alone as 
K = —df^p/p. The function p(0,w) is determined by the 
metric on the curve at u = 0. We will always choose 
the coordinate v to be dimensionless, so that p carries 
dimensions of length. Here, p(u, v)dv is the infinitesimal 
length along curves of constant u and the width, w, is 
the width of the strip after swelling. 

To quadratic order in the strain, the elastic energy is 

m 



We will first consider the general problem of deter- 
mining the shapes of elastic sheets with inhomogeneous 
patterns of swelling. We review the elasticity of such 
swelled sheets in terms of a "preferred" metric assigned 
to the sheet, and discuss the obstructions to minimizing 
both the stretching and bending energies. Finally, we de- 
rive a simple model to describe the shape of the center 
line of a non-Euclidean strip. 



A. Elasticity as differential geometry 

In this section, we review the elasticity of swelled 
sheets, and establish notation for the rest of the paper. 
Consider an elastic sheet with thickness t that has swelled 
by the position-dependent scalar Vt. For thin sheets, the 
energy can be decomposed into two terms, E = Ec + Eb, 
the in-plane stretching e nerg y Ec ex: t, and the bending 
energy Eb oc t^ [10, [l3, 123 ■ This swelling factor mea- 
sures the change of infinitesimal distances on the sheet, 
and can be converted to a "prescribed" metric tensor, 
gij. When the actual metric (tensor) gij equals the pre- 
scribed metric g^, the in-plane stretching energy, Ec 
vanishes. Therefore, Ec is written in terms of a strain 



Ec = 



tY 



tensor, 7ij = gij-gij 



When 17 = 1, 



m an appro- 



priate coordinate system, and this recovers the standard 
formula for the strain of an elastic sheet 23 1 . 

It will prove convenient to consider a coordinate system 
(u, v) in which the metric of the buckled surface takes the 
form 

guudu^ + 2guvdudv + g^dv'^ = du^ + p'^{u, v)dv'^. (1) 

This can be compared to the metric on a flat disk du^ + 



2(1-1/2 

+(1 - |yW''g'^J^Jlkl 



dudvy^ V {g^-'^ij 



(2) 



where Y is the Young's modulus and v the Poisson ratio 
of the sheet. The bending energy is given by the usual 
expression 



Eb =t^- I dudv^ [4H^ - (1 - >')K] , 



(3) 



where k = y/[12(l — i^ )], H — g^^hij/2 is the mean cur- 
vature and K is the Gaussian curvature. These are both 
written in terms of the second fundamental form, hij — 
N • didjT where N is the unit surface normal. Couplings 
between 7^ and hij cannot exist if the swelling is invari- 
ant under a change in surface orientation, N — >■ — N. We 
will also neglect variations of the sheet thickness with 
swelling. 

The total energy E = Ec + Eb, or some variant of it, 
has been considered by a number of authors 0, |Ij, [Ij- 
[l6|. This functional provides a natural generalization of 
the usual elastic energy to non-Euclidean metrics. 



B. Geometrical constraints 

Minimizing this energy is difficult because the metric 
and curvature tensors are not independent when a surface 
is immersed in three dimensions. In fact, they satisfy a 
set of compatibility relations, which we will briefly review 
in this section. To begin, we write two equations for how 



the frame ((9„r,9t,r/p, N) changes |25| . 



C. Non-Euclidean strip energy 
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huv/p 
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-huv/p 



and 



Rn 



duP hvu 

-duP hyy/p 



(4) 



(5) 



(6) 



(7) 



These follow from the definitions of the metric and second 
fundamental forms in our coordinates. 

These two matrix equations are only compatible if 
9ijR„ — dyTlu + [R-MiR-uli where [,] is the commutator. 
This results in three equations relating the second funda- 
mental form to the metric, the Gauss-Codazzi-Mainardi- 
Peterson equations J25i], 



p 

du { -^^ ] = huudup + dv 1 -^ 

p J \ p 



(8) 



The last equation is the theorema egregium of Gauss, and 
will be used to relate the function p to the functions hij . 
Eqs. dl]) and (HJ can also be used as evolution equa- 
tions to construct isometric embeddings for arbitrary p 
provided that hyy ^ (see, for example, [1||| and [18,]). 

Finally, we will find it useful in Sec. |ll] to relate the 
components of the curvature directly to the mean curva- 
ture, H, given by 2H = g^^hij = huu + hyy/ p^. Using 
this relation in the theorema egregium^ we find 



2pH [^ 
P 



pdlp+{hvv/pf 



(9) 



Therefore hyyjp and huv lie on a circle and 
huu = H + viJ2 — X sin a 

\/H'^ — K cosa. 



P- 
h„ 



(10) 



Here, a tells us how the principal curvature axes are 
aligned with respect to the coordinate system. For exam- 
ple, when a = tt/2, the principal curvature H+y'H'^ — K 
points along the tangent vector in the u direction. 



In this paper, we will consider the lowest order terms 
of the stretching and bending energies in an expansion 
in powers of the ribbon width. Our procedure can be 
justified by a careful perturbation expansion of both the 
stretching and bending energies in powers of the strip 
width w, valid as long as w\/—K <^ 1, where K is 
the typical Gaussian curvature associated with the pre- 
scribed metric. This expansion also yields corrections to 
the strain along the center line due to the bending energy 
and will be explored elsewhere |22l |. 

For our purposes, we need the following fact: as the 
ribbon thickness is decreased, its metric more closely 
agrees with the prescribed metric [Ij, [l^- Therefore, 
we will assume our strips are sufficiently thin that the 
strain along the center line vanishes, implying that the 
Gaussian curvature is equal to the prescribed curvature. 
From the numerics in Sec. IIIH the assumption of vanish- 
ing strain along the ribbon center line is roughly borne 
out: as the thickness decreases, the strain along the cen- 
ter of the strip similarly decreases. According to Fig. 
[71 however, the Gaussian curvature mismatch, implying 
strain, remains along the edge. We will see that even the 
condition that K is isometric on the center line is only 
approximately true for the closed strips we are able to 
simulate numerically. Nevertheless, the small remaining 
strain leads to only small quantitive errors. 

The size of this error can be estimated if we note that 
Eqs. ([H can be used to obtain the second fundamen- 
tal form for a fixed metric. For clarity, we will suppress 
all dependence on u unless otherwise indicated until sec- 
tion IIIBl as all quantities in our analytical model will 
be evaluated on the center line at u = 0. If we suppose 
that the metric deviates from the prescribed metric at 
the center line, we can obtain an upper bound on the 
strain. To do so, we suppose a form for the strain for 
which 7„„ = juv = 0. Then jyy ~ p^ ~ p^. Since the 
Gaussian curvature determines d^p, ^yy ^ A7 [u/pY, 
where A7 describes the leading order deformation of the 
strain. Substituting this into the stretching energy in Eq. 
([2]), we obtain a rigorous upper bound on the stretching 
energy 



p^ 



(11) 



where we have dropped a very small numerical factor 
depending on v. This result suggests that the magnitude 
of the energy in the strain can be controlled by adjusting 
the ribbon width or, alternatively, its thickness. 

Similarly, when K — K on the center line, similar con- 
siderations show that Ec ^ Ytw^ jp' ISr{^ where A7 de- 
pends on the third derivative of p at the center line. From 
Gauss' theorem, egregium^ this is linked to derivatives of 
the curvature with respect to u. 

We will approximate the bending energy by its lead- 
ing order power of w, which should be valid when w is 



sufSciently small. For the bending energy, this yield 

t^KW 



Ef 



dv p{0,v)H^{0,v), 



(12) 



where the factor of w accounts for the integral over u. 
Thus, using Gauss' theorem egregium again, this only 
depends on d^p along the center line. 

This suggests that the lowest order terms in the en- 
ergy are minimized when the curvature along the center 
line minimizes Eq. (|12|) subject to the constraint on the 
Gaussian curvature. The higher-order terms 

This type of approximation has also been used to de- 
scribe narrow developable strips. For instance, to ob- 
tain the Sadowsky functional 27], we would minimize 
Eq. p^ subject to the constraint that K = on the 
center line. In Sec. IIIIl we will compare our results to 
a numerical minimization of several closed strips. This 
will both justify our assumptions as well as confirm the 
predictions obtained from our calculations. 

For infinitely thin open strips of sufficient narrowness, 
it is possible to find an isometric immersion such that 
both H = and duH = along the center line f2g|. 
Most of our numerical strips do not even have _ff = on 
the center line; in fact, we will derive conditions in Sec.HIl 
under which it is possible to at least ensure that H — 
on the center line. 



II. CLOSED STRIPS 

We consider strips that are topologically closed and 
unknotted, such as strips resulting from punching a hole 
into the center of a disk or swelling a cylinder. Though 
we will formulate the problem quite generally at first, we 
quickly find ourselves unable to generate sufficient and 
necessary conditions for a strip with a prescribed metric 
to close. To remedy this, we will restrict consideration to 
a large family of strips that yield an analytically tractable 
way to enforce both the prescribed Gaussian curvature 
and the closure constraint. We will therefore generate 
sufficient, but not necessary, conditions for a strip of this 
type to close. 

These strips will have the property that their surface 
normal on the center line always points along a cone. 
Though many of the results in this section are not a pri- 
ori valid for all closed strips, we will show in Sec. IIIIl 
that the numerically minimized strips seem to have this 
property as well. From this assumption, we derive sev- 
eral predictions. In particular, a key result of this section 
is that the strip center line can be minimal (i.e. H — 0) 
if a "resonance" condition 



duP^ + pdlp ~ 1 
2(a„p±l) ' 



(13) 



relating the number of wrinkles m to the derivatives of 
the metric at the center line, is met. We corroborate this 
somewhat surprising result and others numerically in the 
next section. 



A. The class of closed, conical strips 

Closed strips cut from an unswelled, fiat disk have a 
natural coordinate system (r, 0) , where 9 is the azimuthal 
angle around the strip. In the swelled case, it is natural 
to force the coordinate v to be 27r-periodic as well. Once 
we determine p, d^p, &t.P, and hij on the center line at 
u = that result in a closed strip, Eq. (j4]), which governs 
the evolution along lines of constant v, will preserve the 
periodic structure of the strip away from the center line. 

We first note that Eq. ([5]) determining the evolution 
of the strip frame along v can be recast as an equation 
for the components of the two tangents and the surface 
normal in space along the center line. Taking the x- 
component. 




R„ 



' 9„r • X " 

p 
N-x 



(14) 



Similarly, we can replace x with y or z. The skew- 
symmetric matrix R„ is an infinitesimal rotation that 
depends on the second fundamental form along the strip 
center line, which itself is a function of v. Therefore, Eq. 
(HU can be recast as the equation for a particle in a mag- 
netic field. To make this mapping explicit, we identify 
the V coordinate with time and define a "time-dependent 
magnetic field" 



B{v) = - 
and the vector 

Vx(w) = {dur-x)i + 



'-i + huvj - dupk, 



(15) 



P 



j + [N-x]k. (16) 



The index on Vx refers to the fact that it arises from 
taking the a;-component of the frame in Eq. (ITi|) . The 
curves Vy , Vz are defined similarly. The i,j,k introduced 
here label the 9„r, ^^, and N components in the frame 
in Eq. (1141) and are not the directions x,y, z of our three 
dimensional Euclidean space. Eq. (|14p now reads 



9„Vx = B X Vx 



(17) 



So far, the steps we've taken here in no way restrict the 
generality of our discussion. Requiring that our surface 
is a closed strip is equivalent to finding conditions on B 
such that Vi{v) is a closed curve and that /^ p{v) [vi(w) • 

j]dv — for i — x,y, z. Supposing this could be done, the 
problem would be to minimize the energy with respect 
to hij and p on the set of all closed solutions. However, 
this is related to the open Fenchel problem of finding 
conditions on the curvature and torsion of a space curve 
in order for it to close. Instead of finding all possible 
closed solutions, we will construct an ansatz for a class of 
solutions that admit a relatively simple characterization. 
We will restrict ourselves to strips whose normal vec- 
tor, on the center line, is constrained to lie on a cone. For 



convenience, we will refer to these as "conical strips." 
Though our motivation for considering these solutions 
was originally to simplify our analysis, this choice of 
ansatz turned out to be particularly apt, as will be shown 
in Sec. IIIII In appendix|Xl we derive expressions for con- 
ical strips with general mean curvatures along the cen- 
ter line. We summarize the results and notation here. 
First, we define sin^ = H/\/ H'^ — K and recall that 
H = huu + h.uvl (? and K = —d'^pjp on the center line 
from Sec. II Bl We also have the angle a between the 
principal curvature axes of the strip and the axes of the 
uv coordinate system. 

The final result for the tangent vector of the center line 
of the strip is 

-^ = - sin^ I [sin(M/ -I- /3)x -H cos(VF -I- f)y\ 
P ^ 



+ cos^ ^ [sin(H/ - P)x + cos{W - /3)y] 
— sin0cos/3z, 



(18) 



with a constant cone angle (j) between N and z, satisfying 



tan( 



\Jpdlp cos a 



cos/3cos^(9„p- 9„/ 



(19) 



with W{v) = seccj) ^^ dv'[dv' P{v') — dup{v')] + WQ, where 
W{) is an arbitrary constant, and tan/3 = (sin a — 
sin ^)/ cos a = h^vl{phvu)- These expressions can be 
written directly in terms of the functions p and /i^ that 
characterize the shape of the strip on the center line. 

We now derive conditions under which a conical strip 
closes. We first note that the tangent vectors can only 
be periodic if A/3 = [^(27r) - /3(0)]/(27r) and AW = 
[W{2n) - W{0)]/{2tt) are integers. Note that if we fix 
< (j) < tt/2 then the magnitude of the term with 
cos^((/)/2) in Eq. (|18p is larger, so that the number of 
times the strip wraps around will be set by the period 
oi W — P rather than W + /3. Therefore, we interpret 
AM^— A/3 as the number of times the strip wraps around 
the z-axis. In particular, this implies that if we were to 
cut a closed strip from the edge of a disk, we should re- 
quire AW—Af3 = ±1; other choices of APF— A/3 could be 
used to represent knotted strips. In addition, the integer 
A/3 counts the number of times that the tangent vector 
of the center line crosses the z = plane, since it does so 
whenever /3 is an odd integer multiple of 7r/2. Therefore, 
we interpret A/3 as the number of wrinkles of a closed 
strip. Finally, we define Ap = [J^^ dv dup{v)]/{2n). 
From the definition of W, 



AW 



(A/3 - Ap) sec ( 
A/3-Ap 
A/3±l ' 



(20) 



and thus the two choices of AW for a given A/3 mean that 
we have two possible cone angles (j) for a given metric that 
can lead to closed conical strips. This gives us a criterion 
for the existence of such a strip, since Ap and A/3 must 



be chosen so that < cost 
particular, this implies 



< 1 (since < < tt/2). In 



A/3- Ap 



(21) 



The above conditions ensure that dyv/p is periodic. As 
a final constraint, dyV must integrate to zero as v goes 
from to 27r to ensure that r is periodic as well. 



B. Some examples 

It is instructive to consider some specific examples to 
understand the formulas derived in the last section. Con- 
sider first conical strips with A/3 = 0. This class of shapes 
contains all of the axisymmetric strips, for which (3 is con- 
stant in V. With this, < Ap < 1. For axisymmetric 
strips, for which p is also independent of v, this criteria 
is well-known (see the proof in Marder et al. [2g| for an- 
other derivation of this criteria for axisymmetric strips) . 
The class of strips with A/3 = contains, but is not lim- 
ited to, axisymmetric shapes - the function /3 need not 
be constant; all that is required is that /3 be periodic as 
V varies from to 2tt. The condition < Ap < 1 yields a 
limit on how quickly p can grow with u for such a conical 
strip. This has ramifications for a disk with K = — 1/i?^, 
which has a metric with p{r,v) — Rsmh{u/R) |9l. [la|. 
where r is the radial distance from the disk center. If we 
imagine cutting a narrow strip from the edge of the disk 
with center line at r = uq, 9„p = cosh.{uo/R) > 1, which 
violates the constraint for A/3 — 0. 

We next turn our attention to the case of A/3 = 2, 
corresponding to a saddle shape in which ho^w has four 
zeros (two wrinkles) . Then solution in this class of strips 
must have 1 < Ap < 2 or — 1 < Ap < 2. If we specialize 
again to a strip cut from a constant negative Gaussian 
curvature disk so that p(u, v) — i?sinh((M -I- uq)/R), p = 
Rsinh{uo/R) and Ap = cosh{uo/R), and we find that 
a solution of the first type is possible when < uo < 
i?cosh~ (1) or < uo < 1.32i?. At the limiting case, 
cos = 0,1 and the mean curvature along the center line 
diverges, as can be readily seen from Eq. (jASp . 

However, this should not be interpreted as the maxi- 
mum radius of finding a saddle-shaped isometric embed- 
ding. In fact, arbitrarily large isometric embeddings of 
disks with K = —1/R^ exist ^2Q|]. If we were to cut a strip 
from the edge of such an isometric embedding with radius 
larger than i? cosh" ^(2), the strip center line may still 
have two wrinkles (alternatively, four zeros for /ii,„). It 
could not be a conical strip, however. Consequently, it is 
likely that minimal energy strips deviate from the shapes 
we would predict for conical strips, especially when the 
mean curvature is large. 

We now return to conical strips with H = along their 
center line. In that case, using Eq. (|A5p . we have 



duP - dy/S = - cot 4iyJ pdlp 



(22) 



at every v. For simplicity, we will specialize to the case 
that p is constant in v along the center line, though by 
replacing these quantities by their averages over v we 
could lift this specialization. In this case, Ap = duP 
which implies that 0^/3 — A/3. Squaring both sides and 
using cot^ (j) = (^^p^tyi-iAp-Apf derived from Eq. ^, 
we find 



(A/3±lf-(A/3-Apf =pdlp^ 



which implies 



A/3 = 



{duP)^ + pdlp-i 

2(9uP±l) 



(23) 



(24) 



Consider again the disk with constant negative gaussian 
curvature —1/R^. For a strip with center line at radius 
Mo, p — Rsmh{uo/R), duP — coshiuo/R), d^p = p/R^ 
and the resonance condition above is satisfied when A/3 = 
cosh(iio/i?) =F 1- This generalizes the condition found by 
one of us numerically 118]. For the class of solutions 
considered there, this condition could also be derived by 
Fourier analysis. 

Again, it need not be the case that there are no closed 
strips with H = along their center lines when this con- 
dition is not satisfied, only that there are none in the class 
of closed, conical strips. Nevertheless, our results do have 
ramifications on the shape of minimal energy strips: we 
will see in the next section that we can satisfy both the 
bending and stretching energies to lowest order in pow- 
ers of w/R at one of these resonances. Therefore, in the 
regime of thickness we are primarily concerned with in 
this paper, dy/3 = cosh(uo/i?) T 1 is sufficient to demon- 
strate that H — along the center line. The numerical 
evidence we will present in Sec. IIIII suggests that, at the 
very least, embeddings with zero mean curvature on the 
center line do not exist for all choices of metric. 



C. Energy minimization 

We now seek to minimize the bending energy along the 
center line of a closed, conical strip. The energy in Eq. 
p2)l would be minimized when H — along the center 
line. Happily, such a condition can be satisfied within the 
class of closed, conical strips when Eq. ([M)) is satisfied. 
In this section, we are interested in the minimal energy 
shape close to i/ = as well. 

When H = {) along the center line, (3 = mv from Eq. 
([22|) . Near an i/ = resonance, we can look for an 
approximate minimum of the form (3 — mv -\- /, where 
(3 = mv is the solution for the H ~ strip with m 
wrinkles. Using Eq. (IA5[) . we expand H in powers of / 
to find 



H : 



Vpd^P 



tan( 



ta.n(j}{dup-m) \Jp&ip 



dj 



sm{m,v) 



(25) 



This allows us to expand Eq. ([12]) to quadratic order in 
/, from which we find the Eulcr-Lagrange equation 



dy 



dyf 



which has solution 



/^ 



0, 



sin(2777,w) 

2m ^ ' 



(26) 



(27) 



where ^ is a constant independent of v. Looking at Eq. 
(|A5|) for the mean curvature, we see that dv/3 always 
appears in the combination d^p — dy/3. Therefore, we 
can interpret the first term of Eq. (P7|) as a shift in duP, 
which must be accompanied by a w-dependent correction 
proportional to sin(2TOw). Therefore, when duP is tuned 
away from a resonance by an amount A/2, the function 
/3 oscillates with a frequency 2m. 

We can use this near-resonance solution as the basis 
for a general conical strip ansatz valid for all metrics: 
(3 = mv + Asm{2mv)/{2m). For sufficiently small A, 
this ansatz does not modify the location of the zeros of 
sin/3 nor does it alter dy(3 at those zeros. It is straight- 
forward to show that the zeros of hyy coincide with those 
of sin/3. At those zeros, dy/3 = m, — A. For H to remain 
finite at those zeros, inspection of Eq. (IA5|) implies that 



dyf3 — duP + cot (j)^J pd'^p when sin/3 = 0. Therefore, we 
conclude that A — duP + cot (f)^ pd^^p — m. 

Finally, it is necessary to confirm that this ansatz 
yields a closed strip. Since duP is independent oiv,W = 
sec (p[duP v — f3{v) + l3{0)] + Wo. Examining the form of /3, 
both W and /3 change sign when v -^ v + ir/m. Inserting 
this into the equation for the tangent vector and inte- 
grating, we find that Jq^ dv p dyv/p = r(27r) — r(0) = 0. 

To summarize the construction in this section, assume 
that we are given some p{u,v). We would like to con- 
struct an ansatz with m wrinkles. Next, we calculate 
ta,n<j> using Eq. ([^0]) . There are two choices, which are 
perhaps not equally good. Once this choice is made, it is 
straightforward to calculate A, /3, W and then we have 
the expressions for the frame on the strip, including Eq. 
P^ for the tangent vector as well as the other compo- 
nents, whose expressions are in Appendix |^ If we in- 
tegrate the tangent vector, we may generate the space 
curve of the center line. To construct the surface, we 
integrate Eq. ([5]) in a small neighborhood of the center 
line. We depict center lines of strips with m = 3, 4, 5 in 
Fig.H 



III. NUMERICS 
A. Numerical model 

To corroborate our analysis, we have implemented the 
bead-and-spring model of inhomogeneous swelling devel- 
oped by Marder and Papaniolaou [ig|. An unswelled 




to the approximate equation 




FIG. 2. (color online) Analytical predictions for the shapes 
of three-, four- and five- wrinkle strips with H = along the 
center line. The center line is the dark black line, and the 
surface displayed is constructed by extending along the vector 
duT from the center line. Here p — Rsmh{uo/R) and uq, R are 
chosen to satisfy the resonance condition cosh(uo/-R) = m — 1, 
when m is the number of wrinkles. 



sheet is represented by a stacked pair of close-packed tri- 
angular lattices with a periodic length of 200 beads and a 
width of 20 beads (thus, 8000 particles). Nearest neigh- 
bors in the lattice are connected by harmonic springs 
(35000 total bonds) with bond energy between two lat- 
tice sites a and b given by 



^hond — (-. 



(AX)^ 



"■ab 



(28) 



9ij 



(29) 



where |AX| is the bond length. The swelling is incorpo- 
rated by setting the equilibrium bond length according 



where Ax* are the components of the vector connecting 
beads a and h in the unswelled membrane and gtj is the 
prescribed metric. This approximates the integral that 
measures the prescribed distance between points a and 
h. The prescribed metric used in Eq. ([^ is given by the 
average of the continuum prescribed metric evaluated on 
the lattice points a and b. Alternatively, we could eval- 
uate the continuum prescribed metric somewhere along 
the bond; as long as the average bond length is small 
compared to i?, the details of how we evaluate the pre- 
scribed length between neighboring points a and b do not 
change our results. 

As shown in Ref. [16|, Eq. ((28)) is compatible with 
a continuum limit using the strain 7^ = gij — gij. 
Note that the 2D strip energy consisting of a sum of 
stretching and bending energies is recovered from this 
3D model when the thickness is small. In particular, 
we note that the bending energy arises naturally from 
the differential stretching and compression of the two 
layers of the sheet. The energy is minimized using ei- 
ther a conjugate-gradient or Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) algorithm, as implemented the GNU 
Scientific Library (GSL). In general, we find that strips 
develop "creases" often as the minimization procedure 
proceeds unless the strips are sufficiently thick and the 
change in the degree of swelling sufficiently small. As 
a practical matter, therefore, we start our minimization 
with thick strips, t « 0.025u;, and small prescribed cur- 
vatures K K, — 10~^w~^. We then increase the curvature 
in a sequence of steps to the desired curvature (up to 
y Kw ~ 1.3) and, once this sequence is complete, de- 
crease the thickness to the desired thickness (down to 
t w Q.QQAw). If we decrease the thickness first, the strips 
often develop creases, though as long as the strips remain 
smooth during this process, the energy minima we have 
found appear to be quite robust against changes in the 
sequence of steps taken. Finally, it is difficult to main- 
tain a smooth strip shape, even at fixed swelling, once 
the thickness becomes too small. Therefore, we are not 
able to probe the limit of truly thin strips; consequently, 
all our strips seem to have some residual strain, even in 
the ground state. 

We find that strips with different numbers of wrinkles 
are metastable, and have explicitly seen this up to five 
wrinkles. Therefore, to select the number of wrinkles, we 
bias the minimization by choosing an initial strip shape 
with a prescribed number of wrinkles m. Typically, we 
use the graph z(r, 0) = ho cos{m9) as an initial condition 
for the lattice, where r and 9 are polar coordinates in the 
xy plane. We have also tried adding additional Fourier 
modes in the initial conditions and have tried amplitudes 
in the range O.Oli < ho < lOOi for sheets of thickness 
t. In all cases we have tried, the number of wrinkles is 
primarily chosen by the dominant mode in z{r, 9), though 
there is a small bias toward lower energy structures with 
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FIG. 3. (color online) A typical minimal energy strip with 
four wrinkles and prescribed metric pi, with K — — (l/9)«;~ , 
uo = (5/7r)iu and t = O.Olui found by numerical minimization. 
Below, the mean curvature of the center line is displayed. 



two wrinkles, and all strips with the same number of 
wrinkles are nearly identical. If no bias is chosen, the 
minimization will take a long time to break the symmetry 
and will go into a state with two or three wrinkles. 

Finally, since we are interested in verifying the 
assumptions and predictions of our analytical calcu- 
lations here, we will focus on metrics with p{u, v) 
independent of the longitudinal direction v. We 
have explored two prescribed metrics in detail; (1) 
Pi{u,v) = Rsmh{u/R + Uq/R) and (2) p2{u,v) = 
R cosh{u / R) / {2tt) + riRsmh{u/R). Both of these met- 
rics have K = —1/R^ everywhere. The first of these 
corresponds to cutting a strip of width w from the 
edge of a buckled disk at radius uq — ■w/2. The reso- 
nance condition for the center line of this metric gives 
cosh(uo/i?) = 771—1. For 777 = 3, this requires a degree 
of swelling which is difficult to achieve numerically at low 
thicknesses. Therefore, to study the resonance condition 
numerically, we use p2(uo) for various values of rj. From 
the resonance condition in Eq. (j24p . we have four solu- 
tions 



1] = m + (777 ± 1) W 1 



1 



(27r)2(m±l)2 



(30) 



2777 ±1 



1 



2(27r)2(777±l) 



+ 0((m±l)-') (31) 



or,, = 7r7-(r77±l)Jl-^^-^^^-^^ 



(32) 



Tl 



Oiim±ir'), (33) 



2(27r)2(r77±l) 
where the ± must be chosen consistently in each equa- 



tion. With the swelling factors accessible to our numer- 
ics, we are able to probe the branch of resonances with rj 
near -f 1 for 2, 3, and 4 wrinkles. 

Geometric quantities such as H, K for our numerically 
minimized surfaces are computed from an interpolation 
of the points on the mid-surface between the two bead- 
spring sheets. 



B. Results 

As an example of a minimial energy surface from our 
numerical calculations. Fig. [3] shows a four wrinkle strip 
with metric from pi and t — O.Olw with prescribed Gaus- 
sian curvature, K = — (l/9)fi;~2. In Fig. |4l we plot the 
Gaussian curvature along the center line for strips for 
three thicknesses t/w — 0.02, 0.015, and 0.01 and other- 
wise identical parameters. As expexted, along the cen- 
ter line K approaches its prescribed value as the thick- 
ness decreases. It is interesting to note that even for 
the thinnest strips we have tried, the Gaussian curva- 
ture always displays localized regions of stretching. Nev- 
ertheless, our results are apparently robust against this 
localized residual strain. 

In Fig. [5l we plot the three components of the normal 
vector on the center line for the thinnest of the strips 
in Fig. m The x and y components oscillate (with am- 
plitude close to 0.25) around zero while the z compo- 
nent is very nearly constant. To quantify this, we com- 
pute the average A^^ = (N • z) sa 0.950 and the variance 
[{{N ■ z - iV2)2)]~i/2 sa 0.01. The analytical angle coscj) 
is predicted from Eq. (|20)) to be 0.951. 

Similar levels of agreement are obtained for other 
choices of metric, curvature and thicknesses, though thin- 
ner sheets have a closer agreement. We tested 2-, 3-, 
and 4-wrinkle strips of thickness t/w = 0.025 with met- 
ric from pi (with uq/w = b/n as are all strips with pi 
in this paper) starting from R — 150 and increasing 
the curvature down to i? « 20. Fig. [B] shows the re- 
sulting average cos </> on the center line as a function of 
duP — cosh(uo/i?), where the error bars show the vari- 
ance of cos (p. Our theory and Eq. pO| predict a linear 
relationship, with zero variance. This is observed up to 
fairly high curvatures - at the curvatures where cos 4> 
begins to oscillate, we observe an intriguing symmetry- 
breaking from strips with Dmd symmetry to strips with 
Cmv symmetry, and the wrinkles curl into tubes. The 
explanation of this transition is outside the scope of this 
paper, but we display this transition in 2-wrinkle strips 
in the inset of Fig. [S) We note that in this symmetry- 
breaking transition, as in the minimization without a bias 
towards a set number of wrinkles, the system seems to 
take a long time to find a minimum, often oscillating be- 
tween different states. 

In Fig. [3 the distribution of Gaussian curvature is 
plotted for three thicknesses {t/w w 0.008, 0.004 and 
0.0015) and Kw^ = —0.029 for prescribed metric pi. A 
boundary layer is evident in all the images. In both cases 




FIG. 4. Gaussian curvature along the center line for the strip 
shown in Fig. Owith thicknesses t/w = 0.02, t/w = 0.015 and 
t/w = 0.01. As the thickness is lowered, the Gaussian cur- 
vature approaches the prescribed curvature (straight dotted 
line). 
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FIG. 5. The three components of the normal for the strip of 
Fig. [3] along the center line, which makes a nearly constant 
angle with the z-axis. The average z-component is (N • z) — 
Nz ~ 0.950 compared to the predicted value cos(^ = 0.952. 



The variance [((iV ■ z - N^y)] 



-1/2 



0.01. 



the degree of stretching decreases near the center of the 
strip as the thickness decreases. 

We further test our model by comparing the numeri- 
cally obtained H along the center line with our analytical 
calculations. In Fig. [3 we show the specific case with 
prescribed metric pi, t/w = 0.006 and K = —0.029. The 
solid lines show H/ ^J\K\ along the center line for strips 
with different numbers of wrinkles from m = 2 to m — b. 
For the purposes of comparison, we rescale v to contain 
one period of the strip, I-k jm. Though these curves agree 
favorably in magnitude with those predicted by our the- 
ory, there is some discrepancy in the precise strip shape. 
This discrepancy may be caused by the failure of the 
Gaussian curvature to agree with the prescribed curva- 
ture exactly or perhaps by the fact that the analytical 
center line may approximate a curve on the strip that is 
not exactly in the center. 

A more stringent test of the analytical theory is the 
surprising resonance condition of Eq. (j24l) . Since reso- 
nances are at relatively high curvatures it is more difficult 



FIG. 6. (color online) Comparison of the cosine of the angle of 
the normal vector to the z-axis on the center line (points) on 
numerically minimized strips to the prediction of Eq. (|20p 
(solid line) for strips with metric from pi, with thickness 
tjw = 0.025 and parameter mo/w = 5/7r for 2-, 3-, and 4- 
wrinkle strips (red (left), green (middle), and blue (right) 
curves, respectively). The x-axis is d^p = cosh(uo/-R), so 
prescribed curvature increases to the right. Error bars mark 
the variance of cos(/) on the center line as v runs from to 
271. Beneath the plot is a sequence of 2- wrinkle strips show- 
ing the symmetry-breaking transition that causes the conical 
strip ansatz to breaks down at large curvatures. 
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FIG. 7. (color online) Gaussian curvature for strips with 
t/w ^ 0.008, 0.004 and 0.0015 for m = 2 and m = 3. These 
strips have a prescribed metric of pi with Kw — —0.0324. 
Only the middle third of the strip width is shown. Darker 
regions have Gaussian curvature nearer to the prescribed cur- 
vature and so have less strain, while cyan or yellow (lighter) 
regions have more strain. 



to achieve low thicknesses without creases using metric 
Pi(m), we focus on the metric p2{u) with various values of 
■q. We choose t/w = 0.00375 and Kw^ = -0.01; though 
very thin, even at this t/w there is some deviation of K 
along the center line from its prescribed value. In Fig. [HI 
we show the numerically measured variance of mean cur- 
vature and the analytical predictions from the ansatz in 
Sec. Ill CI for such strips. It is clear that the shape of the 
variance as a function of rj is fairly well-mimicked with 
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FIG. 8. The mean curvature HR along the center line of a p\ 
strip with t/w = 0.006 and Kw'^ = -0.029 with 5, 4, 3, and 2 
wrinkles (solid) compared to theoretical predictions (dotted) 
along the center line. 



0.8 


\ 








0.6 


• -•^ 




/ 




0.4 


> \ / 




/ 






* \ / 




>«..„.••" 




0.2 


» //" 










r/ 




Nl/, , , 


1 1 1 1 1 1 1 



1.5 2.0 2.5 3.0 3.5 4.0 
QuP 



vww 




1.010 1.020 

ri 



1.030 



FIG. 9. (color online) Comparison of numerical mean curva- 
ture variance y''{{HR)^) (dashed lines are meant to guide the 
eye between the data points) to analytical predictions (solid) 
for p2, t/w — 0.004 as a function of rj for 2-, 3- and 4- wrinkle 
strips (green, red and blue, respectively) - note the observed 
resonance minima at roughly rj — 1.0120, rj = 1.0060 and 
rj = 1.0042, compared to analytical predictions of 77 = 1.0128, 
77 = 1.0063 and rj = 1.0042, for 2-, 3- and 4-wrinkle strips, 
respectively. 



77 above resonance, but there is less good agreement be- 
low resonance. In fact, the numerically minimized shapes 
tend to look nearly flat at 77 = 1, as opposed to the ansatz 
which predicts a singular strip. For completeness we in- 
clude Fig. (TU] which shows this data for the pi strips 
shown in Fig. [S] - the agreement appears to be worse as 
a result of the thickness being around t/w = 0.025. It 
does appear for 3- wrinkle ribbons that the point at which 
the numerical data begins to differ from the predictions 
occurs near the same values of duP where the symmetry- 
breaking transition causes the conical strip ansatz to fail. 
For 2-wrinkle strips, at 77 w 1.0055 and 77 « 1.0010 
below resonance we observe fairly sharp and reproducible 



FIG. 10. (color online) Comparison of numerical mean curva- 
ture variance •^/{(iifi?)^) (dashed lines are meant to guide the 
eye between the data points) to analytical predictions (solid) 
for pi, t/w = 0.025 as a function of duP = cosh{uo/R) for 2-, 
3- and 4-wrinkle strips (green, red and blue, respectively) - 
note the observed resonance minima at roughly dup ~ 1.193, 
dup = 1.997 and duP = 3.073, compared to analytical pre- 
dictions of duP = 1, 2 and 3, for 2-, 3- and 4-wrinkle strips, 
respectively. 



jumps in the variance of the mean curvature - we do not 
currently have a good explanation of this, as the strip 
shapes do not appear to change dramatically over these 
values of 77. 

In Fig. [m we plot {H^{u)/H^) 

J dv p{u)[H^ {u,v)/H^]/ [271 p{u)] over the surface of 
numerically minimized strips with two, three, and four 
wrinkles above, near, and below resonance. Here Hq is 
the average mean curvature in a certain region of the 
strip. There is a clear trend toward extinguishing the 
mean curvature near the center line of the strip when 77 
is tuned near the resonance condition of Eq. ([23]). 

Eq. ([M)) and its application to p2 strips, Eq. ([5^ . pre- 
dicts rather well the observed location of the resonance. 
We find that the center line of a two wrinkle strip will 
have H — when 77 « 1.01275, a three wrinkle strip will 
have H = when 77 sa 1.00634, and a four wrinkle strip 
will have H = when 77 « 1.00422. 

For example, with two wrinkle strips wc find that the 
mean curvature is slightly smaller with 77 = 1.012 than 
near the predicted resonance of 77 — 1.013. We ascribe 
this small discrepancy to our assumption that the Gaus- 
sian curvature along the center line agrees with that pre- 
scribed by the metric. Indeed, thick membranes do not 
show any evidence of the resonance condition at all, nor 
does the Gaussian curvature reach its prescribed curva- 
ture on the center line. It is also interesting to note that, 
below resonance, the mean curvature near the center line 
of a strip is n out of phase with the mean curvature near 
the strip edges. 
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FIG. 11. (color online) Plots of (H^/Hi) ^ J dv p{u) [H'^{u,v)/ Ho]/[2Tvp{u)] for strips of width w in the region ~w/A < 
u < u)/4, where Ho is the average mean curvature within this region. Tlie prescribed metric is p2{u) and the thickness is 
t/w = 0.00375. (a) {H'^/HI) for 77 = 1.02 (dashed), ry = 1.01 (solid) and r? = 1.00634 (dotted) for m = 2. The resonance 
condition yields 77 « 1.01275. (b) (H'^/Hi) for rj = 1.01 (dashed), rj = 1.00634 (solid) and 77 = 1.0042 (dotted) for m = 3. The 
resonance condition yields rj « 1.00634. (c) {H'^/HI) for 77 = 1.00634 (dashed), 77 = 1.0042 (solid) and r} = 1.003 (dotted) for 
7n = 4. The resonance condition yields 77 « 1.0042. Above is the mean curvature of each strip near resonance with dashed lines 
indicating the limits of the plots below; the mean curvature is clearly extinguished (dark regions) near the center of the strip 
but not near the edges. 
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FIG. 12. (color online) Energy for p = pi as a function of _R = 
Xjy —K for 4-, 3- and 2-wrinkle strips (from top to bottom; 
blue, red and green respectively). 



Finally, we plot the energy for p\ and p2 in Figs. [T2l 
and[T31 In both cases, the energy of to = 2 appears to be 
lowest throughout the range of parameters tested, though 
for a large degree of sweUing the energies become nearly 
degenerate. For Fig. [T31 dips in the energy occur at the 
resonances in 77 as we might expect, since the dominant 
term in the bending energy vanishes near the resonances. 
In Fig. [121 the resonances occur in the regions where R is 
small and the energies are nearly degenerate. At large i?, 
the energy scales as XjB? = \K\, which we expect since 
Eb^ H^ (X l/i?2. 

Note that, at a resonance, the energy is dominated by 
higher order terms that we have neglected in our analysis. 
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FIG. 13. (color online) Energy for p = p2 as a function of rj 
for 4-, 3- and 2-wrinkle strips (from top to bottom; blue, red 
and green respectively). 



It is these additional terms that set the energy scale at the 
resonances. It is these higher order terms that seem to 
determine which shapes have the lowest absolute energy. 



IV. CONCLUSION 

In this paper, we have studied the interplay between 
topology and buckling for elastic strips which have been 
swelled inhomogeneously. We have focused on negative 
Gaussian curvatures with metrics that depend only on 
distance from the strip center line, though nothing pre- 
cludes studying more general metrics. We identify a class 
of closed strips for which we can guarantee a prescribed 
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Gaussian curvature; these "conical" strips have the prop- 
erty that their unit normals point along a cone every- 
where on their center line. Using this family of strips, 
we predict the existence of minimal bending resonances 
for which the mean curvature along the strip center line 
vanishes for a prescribed number of wrinkles. This sur- 
prising result was confirmed in numerical simulations. To 
our knowledge, this phenomenon has not been observed 
experimentally in any system as of yet. 

The main surprise from our numerical tests was the 
fact that conical strips were such a good description of 
the numerical strip shape. It is not clear if this agreement 
can persist to very thin strips, nor is it clear if there are 
other, lower energy, strip shapes which do not exhibit this 
property. Strictly speaking, our analytical calculations 
are only valid within this class of strips and, at best, our 
resonance condition is sufficient, but not necessary, to 
find a strip with H = along the center line. Further 
work quantifying closed strip shapes will be necessary 
to understand this fortuitous agreement. For instance, 
there may exist an expansion of the elastic energy around 
conical solutions which should explain the deviation from 
conical strips at very high curvatures. 

Though the experiments of Klein et al. ^ are for wide 
strips with w^J—K ~ 1, our results do shed some light on 
the frustration that may be leading to the high degree of 
wrinkling observed. Unlike experiments [23|, our numer- 
ical minimization shows no indication that the number 
of wrinkles increases with decreasing thickness. Whether 
this is due to the narrowness of our ribbons compared to 
the experiments is not clear. Indeed, it would be inter- 
esting to probe the degree to which strips with various 
numbers of wrinkles are metastable by perturbing the 
minimized strips and re-minimizing. 

Another interesting point about our analytical predic- 
tions is that only the values at u = of the first two 
derivatives of p with respect to u enter the formulas. Pre- 
liminary tests of minimizing strips with an added u^ term 
to pi indicated that the behavior of the Gaussian curva- 
ture on the center line changed by a bit, but the angle 
of the normal was the same and the mean curvature was 
also close, so we expect that more systematic tests will 
confirm our predictions in their regime of validity. 

It is also interesting to consider the results of the nu- 
merical simulations more closely. One of the generic fea- 
tures of all the numerical minima is the persistence of 
strain in the strip, in particular, in fairly localized re- 
gions of the wrinkle. On the boundaries this is plausible 
as one might expect a boundary layer to appear on quite 
general grounds. On the other hand, even along the cen- 
ter line there is a small degree of strain, and we are unable 
to obtain strips that are sufficiently thin for this strain 
to vanish. Though we do not do so here, this residual 
strain can be studied using a more careful perturbative 
expansion of the energy |22l] . 
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Appendix A: Conical, closed strips 

Recall that we wish to find closed solutions of 
9t,Vx(w) = B(i;) X Vx(w), where 



B(t;) = —i + huvj - dupk. 



(Al) 



In this form, the nontrivial dependence of B on u makes 
it difficult to determine the strip shape or to determine 
generally sufficient and necessary conditions for the strip 
to be closed. Of course, B must be 27r-periodic and 
we shall assume this in what follows. Note furthermore 
that this linear, homogeneous equation may be solved 
by Fourier analysis, though the conditions for closedness 
lead to an infinite set of equations relating the Fourier 
coefficients of p and /ly . 

To motivate our construction, consider the special case 
of B with a constant direction but possibly changing 
magnitude. Closed solutions for Vx to this problem ro- 
tate in a plane perpendicular to B, with a frequency that 
is proportional to the magnitude of B. If we rotate B so 
that it's pointing in the fc-axis, then: 

Vx(w) = |wx|[cos(aj'(;)i + shiiuiv)]] 
^ = \B\ 

The class of solutions that we consider build off of this 
case; we allow the direction of B to move on a cone. 
If we allow a rotation of the ijk coordinate frame by a 
w-dependent matrix, 0(w), then B may be transformed 
into a vector that is always parallel to the fc-axis. We will 
consider the solutions in this class for which the direction 
of B lies along a cone with axis the fc-axis. 

To make a bit more sense of this, remember that the 
fc-axis corresponds to components of the normal vector 
- as it turns out, this restriction will force the normal 
vector to point along a cone oriented along the z-axis in 
space. This is why we will call these conical strips. 

Due to this specialization, we are restricting the al- 
lowed values of the second fundamental form to a sub- 
class of possible strip shapes. Our form for B remains 
general enough, however, that many solutions for closed 
strips still exist in the set of conical strips, including some 
with low bending energies. Furthermore, the center lines 
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of strips we find by numerical minimization have this 
property as weU. 

We proceed by rotating B about the k axis with re- 
spect to an arbitrary angle —j3{v). This yields a modified 
equation 9„Vx' = B x Vx', where 



B' = — ( — ^ cos 13 + huu sin /? ) i 
P 



(A2) 



— —hyu cos/3 H sin/3 j 

\ P J 

+ {dy(3 - dup) k. 

We now choose /3 so that the i component of B' vanishes, 
which simplifies the equation somewhat. In particular, 
we choose 

cos /3 — ~ hyu sin /3 



pH — p\i W^ — Ksm.a cos /3 — —p\/H'^ — K cos a sin /3 
(sin ^ — sin a) cos P — — sin /3 cos a 
sin a — sin^ 



tan/3 



cos Of 



(A3) 



where sin^ = H/y/H'^ — K and we've used Eqs. (HI 
The vector B' will therefore point in a direction in the 
(new) j/c-plane. 

Let (/) be the constant angle between B' and the k axis, 



defined by 










tan(/) = 


Bl. 






= 


-Kz 


, cos /3 + ^^^ sin /3 






9„/? - dyfi 








y/pdlp cos a 



cos /3 cos ^{duP — dyf3) ' 

Note that we've applied the relation K = —d^p/ p. 
By eliminating a, we derive that 



(A4) 



H 



tan^ 



2sin^ 



(A5) 



im\(l){dup-dyl3) 
^pdlP 



\JpdlP 



tan ^{duP - 9^/3) 



This expression for the mean curvature tells us that for 
conical strips, = (flat disk limit) or = 7r/2 (cylinder 
limit) is potentially bad. 

As above, a closed particle trajectory requires that the 
particle velocity, Vx', will lie in the plane perpendicular 
to B'. The instantaneous frequency (in v) of this motion 
in the plane is given by a;(w) = (9«/3 — dyp) sec </). To find 
the complete solution, we return to our original unrotated 
frame by rotating the vector Vx' by — /3 to find 



Vx (w) = [sin cos 7 sin W — cos (\> sin 7] A: + [cos 7 cos /3 cos W + cos cos 7 sin /3 sin W + sin sin 7 sin /3] i 

+ [— cos 7 sin /3 cos M^ + cos ^ cos 7 cos /3 sin 14^ + sin (/) sin 7 cos /3] j , ( A6) 



where 7 is an integration constant having to do with 
motion of the fictitious particle along the B' axis and 
W — Jq dv'uj{v') + Wo- The constant Wq sets the ini- 
tial phase of the motion. Thus W is defined so that 
W{0) = 0, or ly = sec 0/0" dv'{dyl3 - dup). 

Recall that the vector Vx contains the x components 
of the entire frame. But we can also use this solution 
for the components Vy and Vz of the strip frame, each 
of which is described by the solution given in equation 
(jA6p and differing only in the integration constants 7 
and W{Q) = Wq. These constants must be chosen to be 
compatible with the orthogonality of the three vectors 
describing the frame. We choose: (1)7 = and Wq — 0, 
(2) 7 = and Wq = 7r/2, and (3) 7 = -7r/2 to define a 
set of three mutually orthogonal vectors at w = 0. Since 
the evolution acts by rotations, these three orthogonal 
vectors remain orthogonal throughout the evolution. 

If we take the i components of the three solutions we 



get 



duY = sin^ I [cos(W^ + P)x - sin{W + /3)y] 

-t- cos^ I [cos{W - (3)x - sin{W - P)y] 
— sin/3sin(/)z. 



(A7) 



If we take the j components of the three solutions, we 
obtain 

— = - sin^ I [sin(W^ -f- /3)x + cos{W + /3)y] 
P ^ 

+ cos^ ^ [s\i\{W - I3)x + cos{W - I3)y\ 

— cos/3sin(/)z, (A8) 

Finally, the k components yield the normal vector 

N = sin (f) sin Wx + sin </) cos Wy + cos (f>z. (A9) 



14 



[1] D. Drasdo, Phys. Rev. Lett. 84, 4244 (2000). 

[2] H. Liang and L. Mahadevan, Proc. Nat. Acad. Sci. 106, [16 

22049 (2009). 

[3] M. Marder, Found. Physics 33, 1743 (2003). [17; 

[4] E. Sharon, B. Roman and H.L. Swinney, Phys. Rev. E 

75, 046211 (2007). [18 

[5] U. Nath, B.C.W. Crawford, R. Carpenter and E. Coen, [19 

Science 299, 1404 (2003). 
[6] A. Goriely and M. Ben Amar, Phys. Rev. Lett. 94, 

198103 (2005). 
[7] E. Sharon, B. Roman, M. Marder, G.S. Shin and H.L. [20 

Swinney, Nature (London) 419, 579 (2002). 
[8] T. Mora and A. Boudaoud, Eur. Phys. J. E 20, 119 [21 

(2006). 
[9] Y. Klein, E. Efrati, and E. Sharon, Science 315, 1116 [22 
(2007). [23 

[10] L.D. Landau and E.M. Lifschitz, Theory of Elasticity, [24 

3rd Ed. (Pergamon, London, 1959). 
[11] J. Dervaux and M. Ben Amar, Phys. Rev. Lett. 101, [25 

068101 (2008). 

[12] S. Nechaev and R. Voituriez, J. Phys. A 34, 11069 (2001). [26 

[13] B. Audoly and A. Boudaoud, C.R. Mecanique 330, 831 

(2002). [27 

[14] B. Audoly and A. Boudaoud, Phys. Rev. Lett. 91, 086105 

(2003). [28 

[15] M. Marder, E. Sharon, S. Smith, and B. Roman, Euro- 

[29; 



phys. Lett. 62, 498 (2003). 

M. Marder and N. Papaniolaou, J. Stat. Phys. 125, 1065 
(2006). 

E. Efrati, E. Sharon, R. Kupferman, J. Mech. Phys. Sol. 
57, 762 (2009). 

CD. Santangelo, EPL 86, 34003 (2009). 
Q. Han and J-X. Hong, Isometric Embedding of Rieman- 
nian Manifolds in Euclidean Spaces, Mathematical Sur- 
veys and Monographs, vol. 130 (American Mathematical 
Society, Providence, RI , 2006) 

E.G. Poznyak and E.V. Shikin, J. Math. Sci. 74, 1078 
(1995). 

J. Gemmer and S. Venkataramani, e-print 
larXiv: 1 005.4442 

B.G. Chen and CD. Santangelo, (unpublished). 
P.G. Ciarlet, J. Elasticity 78-79, 1 (2005). 
M.P. Do Carmo, Differential Geometry of Curves and 
Surfaces, (Prentice-HaU, Englewood Cliffs, NJ, 1976). 
D.J. Struick, Lectures on Classical Differential Ceometry, 
2nd ed (Addison- Wesley, Reading, Mass, 1961). 
E. Efrati, E. Sharon and R. Kupferman (private commu- 
nication). 

M. Sadowsky, Sitzungsber. Preuss. Akad. Wiss. 22, 412 
(1930). 

M. Marder, R.D. Deegan and E. Sharon, Phys. Today 
60(2), 33 (2007). 
E. Sharon, private communication (2010). 



